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The northernmost limit of distribution of Quercus serrata is on Hokkaido, where fragmented popula¬ 
tions occur from the southwestern foothills of the Hidaka Mountains to the Ishikari Plain with an addi¬ 
tional isolated population on the Oshirna Peninsula. These northern marginal populations are considered 
to be the result of past vegetational shifts, but the genetic relationships are not fully understood. In this 
study, we used 11 nuclear SSR loci to genotype 667 individuals from 11 populations of Q. serrata in Hok¬ 
kaido and six populations from the Tohoku region of Honshu. Total genetic diversity in Hokkaido (H E = 
0.693) did not differ from Tohoku ( H E = 0.669). The population differentiation in Hokkaido (G ’ ST - 0.097) 
was not significantly different from that of Tohoku ( G’ S r = 0.090). The STRUCTURE analysis distin¬ 
guished four regional clusters of Q. serrata; the Pacific side of Tohoku, the Japan Sea side of Tohoku, the 
northernmost region, and the Oshirna Peninsula. The findings indicate two ancestral origins from the 
Pacific and Japan Sea coasts of Tohoku in the Hidaka region. The results from principal coordinate anal¬ 
yses indicated that the isolated population on the Oshirna Peninsula is more related to the populations on 
Hokkaido than to the Tohoku populations. 

Key words: Fagaceae, fragmented distribution, northern distribution range, nuclear microsatellite, 
STRUCTURE analysis 


The distribution of plant species is deter¬ 
mined by several factors, such as physiological 
limitations, ecological niche, and geological his¬ 
tory (Heuertz et al. 2004). Plant populations in 
the center of their distribution generally sustain 
stable conditions in a favorable environment for 
the species, while sites at the margin of their dis¬ 
tribution are exposed to less suitable sites for es¬ 
tablishment and to greater competition with other 
species than at their distribution center (Lawton 
1993, Hampe & Petit 2005). Because of these dis¬ 
advantages, local populations at the margins of 
their range of distribution often occur in small, 
isolated populations (Hampe & Petit 2005, Kita¬ 
mura et al. 2015). The present distribution of 
plant species also provides clues to their past mi¬ 


gration history, which is influenced by various 
factors such as changes in the environment, geo¬ 
logical activity, competition with other species, 
and interaction between biotic and abiotic fac¬ 
tors. During migration, distributional changes 
such as expansion, contraction, fusion and isola¬ 
tion inevitably result in changes in genetic varia¬ 
tion by chance and/or by natural selection. There¬ 
fore, the effects of the past reflect on the genetic 
relationships among current populations 
(McLachlan et al. 2005, Magri et al. 2006). 

Quercus serrata Murray (Fagaceae), is wide¬ 
ly distributed throughout the warm temperate 
and lower cool temperate zones of the Japanese 
Archipelago from Kyushu to Hokkaido and the 
Korean Peninsula, often intermixed with other 
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deciduous broadleaved species of trees. It grows 
between 100 andlOOO m elevation and is better 
adapted to a warmer and drier environment than 
its sister species, Q. crispula (syn. Q. mongolica 
var. crispula) in the northern part of its range in 
northern Tohoku and southern Hokkaido. Quer- 
cus serrata is less frequent in central Hokkaido, 
where it is at the northernmost extent of its natu¬ 
ral distribution. 

Uemura & Takeda (1987) investigated the dis¬ 
tribution, ecological characteristics and climatic 
conditions of 129 cool temperate forest species in 
Hokkaido and classified their distribution pat¬ 
terns into 8 chorological groups represented by 
the major plant species; (A) Fagus crenata group, 
(B) Aesculus turbinata group, (C) Arachniodes 
standishii group, (D) Leptorumohra miqueliana 
group, (E) Cacalia delphiniifolia group, (F) Car- 
pinus laxiflora group (G) Ostrya japonica group, 
and (H) whole area group. In this classification, 
Q. serrata belongs to the Pacific Ocean side dis¬ 
tribution of (E), the Cacalia delphiniifolia group, 
with Neillia incisa (Thunb.) S. H. Oh, Acer cissi- 
folium (Siebold & Zucc.) K. Koch, Callicarpa ja¬ 
ponica Thunb., Spuriopimpinella calycina (Max¬ 
im.) Kitag., and Zanthoxylum piperitum (L.) DC. 
This group has a discontinuous distribution on 
the Oshima Peninsula, Hidaka, and Iburi regions. 
Uemura & Takeda (1987) attributed the distribu¬ 
tion of this group to an environment where there 
is > 50 cm of snow depth for 30 days per year. 

Recent progress in DNA markers combined 
with paleobotanical evidence has contributed to 
the reconstruction of migration routes after the 
last glacial maximum in the northern hemisphere 
(Hewitt 2000, Comes & Kadereit 1998, Taberlet 
et al. 1998, Stewart & Lister 2001, Petit et al. 
2002, Lascoux et al. 2004, Lumaret et al 2005, 
Magni et al. 2005, McLachlan et al. 2005, Magri 
et al. 2006). The phylogeography of Japanese 
oaks has been investigated for both the entire and 
parts of their geographic range (Ubukata 2003, 
Kanno et al. 2004, Okaura et al. 2007, Quang et 
al. 2008, Harada & Ikeda 2011, Ohsawa et al. 
2011). Most of these studies have discussed the 
major species, Q. crispula, while few have fo¬ 
cused on Q. serrata. The studies of Kanno (2004) 


revealed that Q. serrata in Hokkaido was mono- 
phyletic and showed low genetic diversity in nu¬ 
clear simple sequence repeats (nSSRs), but the 
genetic relationships among populations of Q. 
serrata in its northernmost range of distribution 
have not been investigated in detail. 

In general, sympatric species of oaks may 
form natural hybrids that result in introgression 
between pairs of species (Whittemore & Schaal 
1991, Lee et al. 1996, lshida et al. 2003). We not¬ 
ed that the closely related Q. crispula grows with¬ 
in or adjacent to Q. serrata in the northern part of 
its range of distribution. A study by Okaura et al. 
(2007) reported that four sympatric species of 
oaks in the Japanese Archipelago, including two 
cpDNA haplotypes shared among sympatric Q. 
serrata and Q. crispula in Hokkaido, shared iden¬ 
tical cpDNA haplotypes, indicating the possibili¬ 
ty of occasional introgression among the species. 

In this study, we focused on the genetic struc¬ 
ture of Q. serrata in Hokkaido northernmost lim¬ 
it of distribution and compared it with popula¬ 
tions from Tohoku. We discuss the consequences 
of introgression from sympatric Q. crispula to Q. 
serrata and elucidate (i) the genetic diversity and 
differentiation at the northernmost distribution 
limit of Q. serrata and (ii) the genetic relation¬ 
ships of the isolated population on the Oshima 
Peninsula. 

Materials and Methods 

Study sites 

Eleven natural populations of Q. serrata 
across the range of distribution in Hokkaido were 
analyzed; two populations from the northernmost 
limit of its range (1-Takikawa, 2-Mikasa), two 
from the Ishikari Plain (3-Oyafuru, 4-Ohyachi), 
six from the Hidaka region (5-Atsuma, 6-Hobet- 
su, 7-Hirotomi, 8-Shin-ei, 9-Utafue, and 10-Kin- 
eusu), and one from the Oshima Peninsula 
(11-Ohnuma) (Fig. 1, Table 1). Voucher speci¬ 
mens were deposited in the herbarium of the 
Graduate School of Science, Hokkaido Universi¬ 
ty (SAP) (Table 1). In two isolated populations at 
the northern limit of its range, Q. serrata is dom¬ 
inant, but mixed with Q. crispula, Acer japoni- 
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Fig. 1. Geographic distribution of Quercus serrata in Japan and Korea showing the 17 populations analyzed in this study (low¬ 
er right). Numbers correspond to Table 1; Northernmost, Ishikari Plain, Hidaka, Oshima Peninsula, and Tohoku. 


cum, Tilia maximowicziana var. yezoana, and 
Sasa senanensis. Two study sites on the Ishikari 
Plain are also isolated. They are lowland popula¬ 
tions mixed with Q. crispula, Acer japonicum. 
Magnolia obovata, M. kobus var. borealis. Aria 
alnifolia, and Ostrya japonica. The Hidaka re¬ 
gion has the major concentration of Q. serrata on 
Hokkaido (Tohyama & Mochida 1978). The con¬ 
tinuous populations in the Hidaka region are 
mixed forests of Q. crispula, Carpinus laxiflora, 
Acerpictum (syn. A. mono), Sasa nipponica, and 


Cacalia delphiniifolia. One, remote, isolated pop¬ 
ulation of Q. serrata is on the Oshima Peninsula, 
where it is mixed with Q. crispula, Acer amoe- 
num var. matsumurae, and Aria alnifolia', the 
ground cover is Sasa senanensis. 

Six populations from the Tohoku region, 
northern Honshu, were selected as references; the 
Japan Sea side of Tohoku (12-Hirosaki, 13-Tsu- 
ruoka, 14-Tokamachi) and the Pacific side of To¬ 
hoku (15-Iwaizumi, 16-Hayachine, 17-Tazawako) 
(Fig. 1, Table 1). 










Table 1. Location and genetic diversity of the 11 natural Q. serrata Hokkaido populations and six references from northern Honshu analyzed in this study. 
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Table 2. Genetic diversity for 11 SSR loci in Q. serrata. 


Locus name 

TA 

H t 

H s 

Ho 

F IS 

F ST 

r st 

bcqm07 

23 

0.866 

0.840 

0.853 

-0.032 

0.029 

0.024 

bcqm27 

10 

0.567 

0.544 

0.553 

-0.040 

0.036 

0.024 

bcqm31 

10 

0.706 

0.681 

0.698 

-0.033 

0.035 

0.039 

bcqm42 

9 

0.722 

0.714 

0.707 

-0.008 

0.012 

0.034 

bcqm50 

6 

0.626 

0.610 

0.601 

-0.002 

0.031 

0.054 

bcqm74 

12 

0.744 

0.724 

0.701 

0.015 

0.028 

0.028 

bcqm76 

27 

0.874 

0.861 

0.831 

0.020 

0.015 

0.041 

bcqm94 

9 

0.402 

0.394 

0.413 

-0.061 

0.019 

0.031 

bcqm325 

14 

0.751 

0.728 

0.740 

-0.034 

0.032 

0.014 

QpZAG36 

10 

0.779 

0.758 

0.775 

-0.038 

0.034 

0.048 

QpZAG46 

16 

0.748 

0.681 

0.661 

0.011 

0.096 

0.145 

Overall 


0.708 

0.685 

0.685 

-0.018 

0.034 

0.051 


TA: Total number of alleles; H T gene diversity in the total population; // S; average gene diversity within populations; 

H 0 \ observed heterozygosity. F IS fixation index, F SJ and R SJ measures of relative gene differentiation among populations defined 
under the infinite allele model (IAM) and the stepwise mutation model (SMM), respectively. 


Sample collection, DNA extraction, and SSR 
analysis 

Inner bark tissue and leaves were collected 
from 16 to 50 randomly chosen mature trees (di¬ 
ameter at breast height, DBH > 20.0 cm) from 
each population (Table 1). Sample tissues were 
stored below 5°C until DNA extraction. 

Thirty milligrams of sample tissue were fro¬ 
zen with liquid nitrogen and ground to a powder 
with a Multi-beads Shocker (Yasui Kikai Co., 
Ltd, Shiga, Japan). DNA was extracted using the 
DNeasy Plant Mini Kit (Qiagen Inc., Tokyo, Ja¬ 
pan). We employed 11 SSR markers, QpZAG36, 
QpZAG46 (Steinkellner et al. 1997), bcqm07, 
bcqm27, bcqm31, bcqm42, bcqm50, bcqm74, 
bcqm76, bcqm94, and bcqm325 (Mishima et al. 
2006) (Table 2). Polymerase chain reaction (PCR) 
was performed with a Multiplex PCR Kit (Qia¬ 
gen), denatured at 95°C for 15 min, followed by 
35 cycles of denaturation at 94°C for 30 sec, an¬ 
nealing at 57°C for 90 sec, and extension at 72°C 
for 1 min, and a final extension incubation of 30 
min at 60°C. Fragment analyses were performed 
with an ABIPRISM 3130x1 Genetic Analyzer 
(Life Technologies Japan, Tokyo, Japan). Frag¬ 
ment sizes were determined by the GeneScan 
500-LIZ size standard (Life Technologies Japan) 
and an ABI GeneMapper ver. 4.1 (Life Technolo¬ 


gies Japan). Genotyping data have been deposited 
in the TreeGenes Database under accession num¬ 
ber TGDR047. 

Data analysis 

The following genetic diversity parameters 
were calculated in FSTAT ver. 2.93.2 (Goudet 
1995): the number of effective alleles per locus 
( Ne '), allelic richness {AR) based on a minimum 
sample size of 16, observed (H 0 ), the total popula¬ 
tion (H t ), within populations heterozygosity (H s ), 
Nei’s unbiased expected heterozygosity ( H E ), and 
the inbreeding coefficient ( F IS ). The significance 
of deviation of F IS from zero was evaluated on the 
basis of 1,000 randomizations. Genetic differen¬ 
tiation among populations was estimated by F ST , 
and R st . Pairwise F sr was also calculated, and its 
significance was tested on the basis of 1,000 ran¬ 
domizations, after Bonferroni correction. 

We used GenoDive ver. 2.0b23 (Meirmans & 
Van Tienderen 2004) to calculate G’ S t (Hedrick 
2005) and compare differences in genetic diver¬ 
sity parameters (Goudet 1995) between Tohoku 
and Hokkaido regions; the test was performed 
under 999 permutations by randomizing the pop¬ 
ulations over the regions. 

The program Arlequin ver. 3.5 (Excoffier et 
al. 2005) was used to test the linkage disequilib- 
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rium (LD) within each pair of loci and popula- 
tions by 10,000 permutations. 

GenAlEx ver. 6.5 (Peakall & Smouse 2006) 
was used to calculate the number of private al¬ 
leles (PA) for the 17 populations and isolation by 
distance (1BD) by simple regression analysis with 
a matrix of pairwise population differentiation in 
terms of F ST l (1 - F ST ). The significance was eval¬ 
uated by the Mantel test based on 999 permuta¬ 
tions (Mantel 1967). Principal coordinate analy¬ 
sis (PCoA) of 17 populations was conducted 
based on pairwise F ST among populations using 
GenAlEx. 

To evaluate recent bottlenecks, we employed 
a two-tailed Wilcoxon signed-rank test using 
Bottleneck ver. 1.2.02 (Piry et al. 1999) under the 
two-phase model (TPM), the infinite allele model 
(IAM) or the stepwise mutation model (SMM). 

STRUCTURE ver. 2.3.4 (Pritchard et al. 
2000) was used to understand the components of 
genetic ancestry among populations. Analyses 
were based on an admixture model and correlat¬ 
ed allele frequencies model, with prior informa¬ 
tion on location (hereafter, the F-model). Simula¬ 
tions were carried out 10 times for each K = 1-10 
value with 70,000 burn-ins and 30,000 MCMC 
iterations. The log likelihood, In P (X\K), and A K 
were calculated in STRUCTURE Harvester (Earl 
& vonHoldt 2012) to estimate the most likely 
number of clusters ( K ). The STRUCTURE re¬ 
sults were summarized in CLUMPP Windows 
1.1.2 software (Jakobsson & Rosenberg 2007). 
Mean F values (equivalent to F ST ) for each cluster 
and the genetic distance among clusters was cal¬ 
culated from the results of 10 runs, considering 
label switching among the runs. 

Results 

Genetic diversity of northern populations of Q. 
serrata 

The genetic diversity of the 11 SSR loci is 
summarized in Table 2. Gene diversity of the to¬ 
tal population (If,) and the average gene diversity 
within populations ( H s ) was 0.708 and 0.685, re¬ 
spectively. The observed heterozygosity (H 0 ) 
ranged from 0.413 at bcqm94 to 0.853 at bcqm07. 


The value of F IS ranged from -0.061 (bcqm94) to 
0.020 (bcqm76) (Table 2). 

Table 1 shows the genetic variation in the 17 
populations. The number of effective alleles ( Ne ) 
ranged from 3.442 (6-Hobetsu) to 3.996 (9-Uta- 
fue) in Hokkaido and from 3.240 (15-Iwaizumi) 
to 4.210 (13-Tsuruoka) in Tohoku. The number of 
private alleles (PA) for the population of 9-Utafue 
was 5, which was the highest in the populations 
we examined from Hokkaido. The allelic rich¬ 
ness (AR) ranged from 5.436 (11-Ohnuma) to 
6.286 (8-Shin-ei) in Hokkaido, and from 5.854 
(15-Iwaizumi) to 6.742 (13-Tsuruoka) in Tohoku. 
H 0 ranged from 0.658 (2-Mikasa) to 0.736 
(8-Shin-ei) in Hokkaido and from 0.617 (16-Hay- 
achine) to 0.699 (14-Tokamachi) in Tohoku; the 
expected heterozygosity (H E ) ranged from 0.652 
(6-Hobetsu) to 0.716 (3-Oyafuru) in Hokkaido 
and from 0.628 (15-Iwaizumi) to 0.721 (13-Tsu¬ 
ruoka) in Tohoku. The F, s values ranged from 
-0.063 (8-Shin-ei) to 0.050 (2-Mikasa). Three 
populations deviated significantly from zero 
(2-Mikasa, 6-Hirotomi, and 11-Ohnuma). 

The proportion of significant (P < 0.05) link¬ 
age disequilibrium (LD) among 55 pairs of loci 
within populations ranged from 0.073 (12-Hiro- 
saki and 14-Tokamachi) to 0.436 (3-Oyafuru and 
4-Ohyachi) (Table 1). Mishima et al. (2006) re¬ 
ported significant LD for 4 pairs of 6 loci which 
we used in this study. However, the observed bias 
of LD was not severe, even if the SSRs were par¬ 
tially linked, as the following results of the 
STRUCTURE analysis was similar to the PCoA, 
which was free from genetic assumptions. 

Genetic differentiation among populations 

The genetic differentiation among popula¬ 
tions, F st , Rst, and G’ ST , was 0.034, 0.051, and 
0.106, respectively (Table 2). The genetic differ¬ 
entiation G’ sr in Hokkaido (0.097) was compara¬ 
ble (P = 0.739) to that of the Tohoku region (G’ ST 
= 0.090). Pairwise F ST values for all population 
pairs, except for 8-Shin-ei to 9-Utafue and 12-Hi- 
rosaki to 14-Tokamachi, showed significant dif¬ 
ferentiation (Table 3). 

The Bottleneck analysis results did not show 
any significant excess of heterozygosity (P < 
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Table 3. Pairwise F ST (below diagonal) and its significance (upper diagonal) among 17 Q. serrata 



1 

Takikawa 

2 

Mikasa 

3 

Oyafuru 

4 

Ohyachi 

5 

Atsuma 

6 

Hobetsu 

7 

Hirotomi 

8 

Shin-ei 

9 

Utafue 

10 

Kine-usu 

11 

Ohnuma 

12 

Hirosaki 

13 

Tsuruoka 

14 15 

Tokamachi Iwaizumi 

16 17 

Hayachine Tazawako 

1 Takikawa 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

2 Mikasa 

0.016 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

3 Oyafuru 

0.033 

0.018 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

4 Ohyachi 

0.032 

0.021 

0.018 

- 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

5 Atsuma 

0.031 

0.020 

0.017 

0.021 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

6 Hobetsu 

0.020 

0.020 

0.025 

0.030 

0.020 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

7 Hirotomi 

0.039 

0.037 

0.026 

0.048 

0.029 

0.032 


** 

** 

** 

** 

** 

** 

** 

** 

** 

** 

8 Shin-ei 

0.045 

0.036 

0.022 

0.042 

0.017 

0.030 

0.016 


n.s. 

** 

** 

** 

* 

** 

** 

** 

** 

9 Utaftie 

0.048 

0.037 

0.022 

0.037 

0.021 

0.027 

0.035 

0.006 


** 

** 

** 

** 

** 

** 

** 

** 

lOKine-usu 

0.033 

0.025 

0.017 

0.024 

0.012 

0.023 

0.027 

0.022 

0.023 


** 

** 

** 

** 

** 

** 

** 

11 Ohnuma 

0.064 

0.050 

0.033 

0.052 

0.043 

0.059 

0.044 

0.040 

0.043 

0.031 


** 

** 

** 

** 

** 

** 

12 Hirosaki 

0.061 

0.051 

0.047 

0.060 

0.046 

0.048 

0.012 

0.024 

0.045 

0.054 

0.054 


** 

n. s. 

** 

** 

** 

13 Tsuruoka 

0.044 

0.037 

0.029 

0.040 

0.030 

0.033 

0.014 

0.011 

0.029 

0.030 

0.041 

0.013 


* 

** 

** 

** 

14 Tokamachi 

0.053 

0.048 

0.041 

0.048 

0.043 

0.039 

0.018 

0.025 

0.043 

0.047 

0.045 

-0.001 

0.017 


** 

** 

** 

15 Iwaizumi 

0.042 

0.048 

0.037 

0.040 

0.038 

0.026 

0.050 

0.037 

0.041 

0.056 

0.058 

0.052 

0.038 

0.043 


** 

** 

16Hayachine 

0.033 

0.036 

0.025 

0.038 

0.039 

0.025 

0.044 

0.031 

0.032 

0.034 

0.049 

0.048 

0.028 

0.047 

0.028 


** 

17Tazawako 

0.037 

0.028 

0.025 

0.038 

0.018 

0.025 

0.037 

0.022 

0.021 

0.028 

0.046 

0.047 

0.020 

0.046 

0.028 

0.018 





Distance (km) 


Fig. 2. Relationship between pairwise differentiation, F ST / (1 — F ST ), and geographic distance (km) among populations of Q. 
serrata ; a) all 17 populations, b) 11 populations in Hokkaido. 
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c) 



Fig. 3. Structure results for 17 populations of Q. serrata. (a) Bar plots for K = 2 to 7. (b) Neighbor-joining tree for each cluster 
when K= 4 along with F values for each cluster, (c) Log likelihood for each K value. Open circles, In P (XLK)\ filled circles, 
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Fig. 4. Principal coordinate analysis (PCoA) results for 17 populations of Q. serrata. Dotted line divides populations between 
Tohoku and Hokkaido regions. 


0.05) under TPM and SMM (Table 1). On the one 
hand, significant P values were detected under 
IAM in populations in Hokkaido, except for 
1-Takikawa, 6-Hobetsu, and 7-Hirotomi. Howev¬ 
er, the IAM may include a type I error (Cornuet 
& Luikart 1996) since its standard deviations 
were greater than SMM or TPM (not shown). 

The 1BD among all 17 populations and the 11 
populations on Hokkaido were significant (Fig. 2, 
R 2 - 0.311 and 0.652, respectively, P < 0.001). 

The STRUCTURE results are shown in Fig. 
3. The rise in the In P (X\K) value slowed when K 
> 4 and reached a plateau at K = 7. The highest 
A K was observed when K = 4, indicating that the 


most likely number of clusters was 4 (Fig. 3c). 
Bar plots of K = 2-7 are shown in Fig. 3a. The K 
= 2 bar plot distinguished the Japan Sea (yellow 
in Fig. 3a) and Pacific sides of Tohoku (blue in 
Fig. 3a), and the Hidaka populations shared both 
clusters. At K = 3, the additional cluster (red in 
Fig. 3a) was assigned to 11-Ohnuma on the Os- 
hima Peninsula. At K= 4, the most likely number 
of clusters, designated the northernmost cluster 
(green in Fig. 3a). Two populations in Hidaka, 
7-Hirotomi and 8-Shin-ei, shared the same ances¬ 
tral cluster with populations on the Japan Sea side 
of Tohoku (Fig. 3a). Furthermore, additional clus¬ 
ters with K = 5 were assignable to the lshikari 
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Plain, namely, 3-Oyafuru and 4-Ohyachi, and to 
the Pacific side of Honshu. The Ishikari Plain was 
distinguished by K = 6. Furthermore, a mutual 
cluster was detected among the northernmost 
Ishikari Plain, Hidaka and Tohoku populations at 
K = l. Consideration should be given to the pos¬ 
sibility of artificial transportation of seeds of Q. 
serrata to Hokkaido during human migrations, 
although we do not currently have evidence to 
discuss the role of humans. 

The F-values for each cluster, analogous to 
the traditional F sr values estimated between the 
common ancestral population and each cluster 
assumed in the F-model, at K = 4 were 0.024, 
0.030, 0.072, and 0.093 (Fig. 3b). The cluster rep¬ 
resenting the Pacific side of Tohoku (blue) had the 
lowest F-value, while the cluster with the second 
lowest F (yellow) consisted of populations on the 
Japan Sea side of Tohoku. The highest F (red) 
was assigned to the Oshima Peninsula. The sec¬ 
ond highest (green) was assigned to the northern¬ 
most region. 

Figure 4 shows the PCoA result by pairwise 
F sr . The first coordinate (COl, explaining 28.61% 
of variation) reflects the geographic locations 
from the northernmost, Ishikari Plain, Japan Sea 
side of Tohoku and Hidaka, to the Pacific side of 
Tohoku. The second coordinate (C02, explaining 
16.90% of variation) distinguishes the Hokkaido 
populations, except for 1-Takikawa and 6-Hobet- 
su, as negative and the Tohoku populations as 
positive values. 

Discussion 

Possibility of interspecific hybridization involv¬ 
ing Q. serrata at the northern limit of its distribu¬ 
tion 

Sympatric species of oaks may show intro- 
gression among species (Whittemore & Schaal 
1991, Lee et al. 1996, Ishida et al. 2003). Ortego 
et al. (2015) studied populations of six species of 
sympatric oaks using a demographic approach to 
determine genetic admixture and hybrid swarms. 
Even though a pervasive haplotype was shared 
among sympatric species, the species maintained 
their morphological and nuclear genetic identity. 


Similarly, the phenotypes of the species in our 
study sites appeared to be stable, since we did not 
find morphologically intermediate individuals. 

The res ults of the STRUCTURE revealed that 
Q. serrata in northern Japan showed genetic dif¬ 
ferentiation among populations (Fig. 3). In con¬ 
trast, Q. crispula in Hokkaido and Tohoku 
showed low levels of population differentiation in 
nuclear SSRs (Ohsawa et al. 2011). That is, popu¬ 
lation differentiation in Q. serrata was not the 
same as in Q. crispula in the northern portions of 
their ranges. Our preliminary analysis using 9 
EST-SSR loci detected a substantial distinction 
between Q. crispula and Q. serrata on Hokkaido 
(Matsumoto, unpublished data); STRUCTURE 
analysis clearly distinguished two species by K = 
2, in which 8% of the individuals of Q. serrata 
and 1% of the individuals of Q. crispula showed 
evidence of natural hybridization with q < 0.9. 
Okaura et al. (2007) shows that Q. serrata in Ho- 
betsu corresponds with 6-Hobetsu in our study 
and Q. crispula in Ohnuma was fixed for the cp- 
DNA haplotype II. Quercus serrata in Ohnuma, 
corresponding to 11-Ohnuma, had two haplo- 
types, I and II, indicating that the composition of 
the haplotypes of sympatric Q. crispula and Q. 
serrata in Ohnuma is not entirely consistent. 

The apparent species delimitation in nuclear 
genomes and the compound structure of organ¬ 
elle genomes among species may be explainable 
by the higher conspecific nuclear gene flow than 
organelle gene flow (Currat et al. 2008, Petit & 
Excoffier 2009). Higher conspecific gene flow in 
our study sites may be supported by phenological 
observations suggesting infrequent interspecific 
crossing, since Q. serrata flowered two weeks 
earlier than sympatric Q. crispula every year 
(Namikawa, pers. obs.). We therefore disregard 
the influence of hybridization with Q. crispula 
and interspecific gene flow between nuclear 
SSRs. 

Another issue is whether the migration and 
expansion of the distribution of Q. serrata was 
independent from Q. crispula at the northern lim¬ 
it of their distribution. From paleobotanical evi¬ 
dence, fossil pollen diagrams showed a rapid ex¬ 
pansion of the range of Quercus in Hokkaido in 
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lowland forests after 8,000 years BP (Yasuda & 
Miyoshi 1998). Morphological similarity of fossil 
pollen of Quercus, however, makes it impossible 
to distinguish between Q. crispula and Q. serrata 
(Hayashi et al. 2010). Although the range of dis¬ 
tribution of Q. serrata and Q. crispula currently 
overlap, their microhabitats differ. Quercus ser¬ 
rata dominates steep, south to southwest facing 
drier slopes where Q. crispula dominated forests 
are absent (Namikawa, unpublished data). We 
therefore do not assume competition between 
these two species, nor any effects resulting from 
their past range expansion. 

Genetic diversity of Q. serrata in the northern 
range of distribution 

Plant populations on the margins of their 
range of distribution generally tend to have low 
genetic diversity (Hampe & Petit 2005, Vernesi et 
al. 2012, Iwasaki et al. 2012, Kitamura et al. 
2015). The southern richness to northern purity 
of genetic diversity in the northern hemisphere 
under the influence of the Quaternary ice age is 
evident (Hewitt 2000). For example, northern 
populations show less genetic diversity in Betula 
maximowicziana (Tsuda & Ide 2005) and Fagus 
sylvatica (Comps et al. 2001) than in their more 
southern populations. A thorough investigation 
of Fagus crenata showed that genetic diversity 
declined northward (Kitamura et al. 2015). Our 
findings, however, revealed average heterozygos¬ 
ity of Q. serrata on Hokkaido (//,, = 0.694) to be 
comparable (P = 0.051) to that of the Tohoku pop¬ 
ulations ( H e = 0.667) (Table 1). Furthermore, the 
northernmost populations (1-Takikawa and 2-Mi- 
kasa) did not show a decline in genetic diversity 
(P = 0.838). Allowing that our study did not cover 
the entire range of Q. serrata, but was confined to 
the northern limit of its distribution, our findings 
are inconsistent with previous studies of major 
deciduous forest trees in the Japanese Archipela¬ 
go, whose genetic diversity shows a decline in 
northern populations (e. g. Tomaru 2013). 

That Q. serrata in Hokkaido does not show a 
northward decline in genetic diversity may indi¬ 
cate that the species is not only the result of recent 
northward expansion, but may have been derived 


from relicts from earlier distributions, which pro¬ 
vided relevant genetic diversity. Recent studies 
suggest the existence of cryptic refugia (Provan 
& Bennett 2008) within northern ranges of distri¬ 
bution (Bhagwat & Willis 2008). The existence 
of relict populations of temperate tree species, 
such as B. maximowicziana (Tsuda & Ide 2010) 
and Fraxinus mandshurica var. japonica (Hu et 
al. 2010), in the Japanese Archipelago, has also 
been discussed. Ohsawa et al. (2011) discussed 
the possibility of cryptic refugia for Q. crispula, 
which is closely related to Q. serrata, in Hokkai¬ 
do. Hence, the existence of northern relict popu¬ 
lations of Q. serrata may explain the ample 
amount of genetic diversity among present popu¬ 
lations. 

The present distribution of Q. serrata is char¬ 
acterized by climatic condition on the Pacific 
Ocean side of Japan. The discontinuous distribu¬ 
tion in the Hidaka, Iburi and Oshima regions 
(Uemura & Takeda 1987) could be due to frag¬ 
mentation of a formerly continuous distribution 
by climatic oscillations (Yasuda & Miyoshi 
1998). Three extremely small and isolated popu¬ 
lations are on Hokkaido; 3-Oyafuru, 6-Hobetsu, 
and 11-Ohnuma. 3-Oyafuru is a small remnant 
population with a small number of juveniles of Q. 
serrata on the low Ishikari Plain. 6-Hobetsu has 
only 30 mature trees and no saplings or seedlings. 
11-Ohnuma is also small and isolated from the 
other populations by more than 100 km (Fig. 1). 
IBD (Fig. 2) indicates little practical gene flow. 
These populations have been more prone to re¬ 
cent disturbances caused by human settlement 
since Q. serrata grows at lower elevations and is 
therefore adjacent to or overlaps with human set¬ 
tlements. The low genetic diversity might be ex¬ 
pected due to the extremely small population size 
and isolation (Young et al. 1996, Lowe et al. 
2005). However, genetic diversity, such as FL E and 
AR, was not significantly different from other 
populations (P > 0.05). We assume that even 
though these populations had experienced popu¬ 
lation reduction and isolation, an ample amount 
of genetic diversity might be secured if the popu¬ 
lation is kept intact after reduction and isolation 
(Kitamura et al. 2015). In addition, Isagi (2004) 
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pointed out that long-lived trees do not show re¬ 
duction in genetic diversity due to recent frag¬ 
mentation because the effects of the fragmenta¬ 
tion on genetic diversity become evident only af¬ 
ter a long period of generational turnover. Trees 
of Q. serrata in the Hidaka region are 55 to 100 
years old (Ishikawa et al. 2015). The maximum 
age for the species may be comparable with other 
species of Quercus, which ranges from 300 to 
400 years (Di Filippo et al. 2015). It is probable 
that the small, isolated populations of Q. serrata 
are the result of fragmentation too recent to be re¬ 
flected in the loss of genetic diversity. 

Genetic structure of Q. serrata in the northern 
range of distribution 

Our study revealed several genetic implica¬ 
tions for the northern populations of Q. serrata. 
Firstly, we recognized ancestral clusters by 
STRUCTURE analysis reflecting geographical 
locations (Fig. 3). The most likely numbers of 
cluster, K= 4, revealed that the lowest (0.024) and 
the second lowest (0.030) F-values represented 
the Pacific and Japan Sea sides of Tohoku, re¬ 
spectively, indicating that Q. serrata in Tohoku is 
ancestral to the more northerly populations (Fig. 
3b). Those clusters were also allocated to Hidaka 
and the Ishikari Plains, although the lshikari 
Plain was assigned to clusters with K > 5. The re¬ 
sults indicate that Q. serrata in Hidaka has con¬ 
served the ancestral lineage related to the Tohoku 
region. Honjo et al. (2004) found that Primula 
sieboldii in the Hidaka region shared a common 
cpDNA haplotype with the Tohoku region and 
corroborated the Hidaka area as being a refugium 
for many temperate species during the last glacial 
period in Hokkaido (Uemura & Takeda 1987). 

In contrast, the population on the Oshima 
Peninsula is highly differentiated for K > 3 (Fig. 
3). Similarly, localized genetic variation in cpD¬ 
NA was found on the Oshima Peninsula (Okaura 
et al. 2007). Euonymus oxyphyllus also exhibits 
local genetic variation on the Oshima Peninsula 
(Iwasaki et al. 2012). It is evident that the differ¬ 
entiation of 11-Ohnuma population can be ex¬ 
plained by 1BD (Fig. 2). However, STRUCTURE 
analysis showed that the cluster representing the 


Oshima Peninsula, red in Fig. 3a, was also pres¬ 
ent in the Hidaka populations. More significantly, 
the PCoA results support differentiation of the 
Hokkaido and Tohoku populations, and the popu¬ 
lation in Oshima Peninsula being classified 
among Hokkaido (Fig 4). Our findings indicate 
that the 11-Ohnuma population is closer to other 
populations on Hokkaido than to those in the To¬ 
hoku region. 
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